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^ i Abstract 

^ ■ The solution of the equations for the generating functions of mul- 

. tiplicity distributions of quark and gluon jets in higher order 3NL0 

\ perturbative QCD is obtained. The non-perturbative effects are dis- 
cussed as well. The results are compared to recent experimental data 

■ on mean multiplicities, their ratio and its slope, higher moments of 

i the distributions. 

in 



o 

oo 



O 
O 



Here we briefly describe recent advances in theoretical understanding of 



O . multiplicity distributions of quark and gluon jets. The extended survey with 



more detailed comparison to experimental data can be found in the recent 
review paper 

^ The progress in experimental studies of properties of quark and gluon 

p • jets is very impressive. Therefore the study of the energy evolution of such 

D , parameters of multiplicity distributions of jets as their average multiplicities 

^ ' and widths becomes possible. It is well known that the average multiplicities 

of quark and gluon jets increase quite fast with energy but their ratio has a 
^ ■ much slower dependence. 

■ The perturbative QCD provides quite deflnite predictions which can be 

confronted to experiment. In brief, the results can be summarized by saying 
that the energy dependence of the mean jet multiplicity can be perfectly fltted 
but the ratio of gluon to quark jet multiplicities can be described within the 
precision of 15-20% only. Moreover, one can understand why next-to-leading 
approximation is good enough for describing the energy dependence, but it is 
not quite satisfactory yet for the ratio value. We show this by presenting the 
analytical expressions. For the corresponding Figures, we refer the reader to 
the review paper 

The theoretical asymptotical value of the ratio of average multiplicities 
equal 2.25 is much higher than its experimental values, which are in the 
range from 1.05 at comparatively low energies of T resonance to 1.5 at Z° 
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resonance. The next-to leading order (NLO) corrections reduce this ratio 
from its asymptotical value by about 10% at energy. The NNLO and 
3NL0 terms diminish it further and show the tendency to approximate the 
data with better accuracy. The computer solution of QCD equations for 
the generating functions has shown even better agreement with experiment 
not only on this ratio but on higher moments of multiplicity distributions as 
well. Being perfect at energy, the agreement in the ratio is not as good 
at lower energies where the theoretical curve is still about 15-20% above the 
experimental one. In other words, the theoretically predicted slope of the 
ratio of multiplicities in gluon and quark jets is noticeably smaller than its 
experimental value. Nevertheless, one can speak about the steady conver- 
gence of theory and experiment with subsequent improvements being done. 
Moreover, it is even surprising that any agreement is achieved in view of the 
expansion parameter being extremely large (about 0.5) at present energies. 

The importance of studying the slopes stems from the fact that some of 
them are extremely sensitive (while others are not) to higher order perturba- 
tive corrections and to non-perturbative terms in the available energy region. 
Thus they provide us with a good chance to learn more about the structure 
of the perturbation series from experiment. 

In the perturbative QCD, the general approach to studying the multiplic- 
ity distributions is formulated in the framework of equations for generating 
functions. Therefrom, one can get equations for average multiplicities and, in 
general, for any moment of the multiplicity distributions p|, |^ . In particular, 
two equations for average multiplicities of gluon and quark jets are written 
as 



Herefrom one can learn about the energy evolution of the ratio of multiplic- 
ities r and of the QCD anomalous dimension 7 (the slope of the logarithm 
of average multiplicity in a gluon jet) defined as 



{nciy))' = J dx^l[K^{x){{nG{y + Inx)) + {nc^y + ln(l - x)) - {ndy))) 

+nfK^{x){{nF{y + \nx)) + {np^y + ln(l - x)) - {nciy)))], (1) 

My))' = J dx-f^K^{x){{nG{y + In x)) + {nF{y + \n{l-x))-{nF{y))). (2) 
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r = 
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Here, prime denotes the derivative over the evolution parameter y = hi{pQ/Qo) 
p, are the momentum and the initial angular spread of the jet, related to 
the parton virtuality Q = pQ/2, (5o=const, i^'s are the well known split- 
ting functions, (n^) and {np) are the average multiplicities in gluon and 
quark jets, {ug) is the slope of {no), nf is the number of active flavours. 
The perturbative expansion of 7 and r is written as 

7 = 7o(l - ai7o - a27o " «37o) + 0(70), (4) 

r = ro(l - ri7o - r27o - ^37o) + 0(7o), (5) 



where 70 = y2Ncas/TT, as is the strong coupling constant, 
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+ 0{y~'), (6) 



Po = {IIN^ - 2n/)/3, Pi = (SIAT, - 19nj)/3,ro = NjCp, and in QCD 
Nc = 3 is the number of colours, Cp = 4/3. 

The limits of integration in eqs. (|1]), (0) used to be chosen equal either 
to and 1 or to and 1 — e~^. This difference, being negligibly small at 
high energies y, is quite important at low energies. Moreover, it is of physics 
significance. With limits equal to and 1 — e~^, the partonic cascade 
terminates at the perturbative level Qo as is seen from the arguments of 
multiplicities in the integrals. With limits equal to and 1, one extends the 
cascade into the non-perturbative region with low virtualities Qi ~ xpQ/2 
and Q2 ~ (1 — x)pQ/2 less than Qo/2. This region contributes terms of 
the order of e~^, power-suppressed in energy. It is not clear whether the 
equations and LPHD hypothesis are valid down to some Qo only or the non- 
perturbative region can be included as well. 

Nevertheless, the purely perturbative expansion (^), with constant 
coefficients Oj, rj and energy-dependent 70 is at work just in the case of limits 
and 1. The values of Oj, rj for different number of active flavors nf are tab- 
ulated in m . At Z'^-energy the subsequent terms in (|^) diminish the value of 
r compared with its asymptotics ro = 2.25 approximately by 10%, 13%, 1% 
for n/ = 4 getting closer to experiment. However the theoretical value of r 
still exceeds its experimental values by 15-20%. At lower energies the value 
of r diminishes due to the running property of the coupling strength. Thus 
these calculations are more sensitive to it than the fractality studies where 
it is almost unnoticed at large bins (e.g., see [§], 0). 
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The energy dependence of mean multiplicities can be obtained Q from 
the definition by inserting there the value of 7 and integrating over 
y. Keeping the terms as small as y"^ at large y in the exponent, one gets 
[|] the following expressions for energy dependence of multiplicities of gluon 
(G) and quark (F) jets 



{no) = Ky"^^""^ e^V pC^^ + 5g(z/)] , 



(7) 



with K an overall normalization constant, C = JANc//3o, and 



C 



2a,C' + %[\n{2y) + 2] 



y 



Po 



with 



W = -2/-"^^%xp [2C^ + 5p{y)] 
f 



C 

Spiy) = Sciy) + — ri + — (r2 + 



(9) 



(10) 



It happens that 2NL0 and 3NL0 terms (contributing to ?/~^^^ and y~^ terms 
in the exponent) are almost constant at present energies and do not change 
the energy dependence prescribed in NLO approximation. It explains why 
the energy dependence is well fitted by both NLO and 3NL0 formulas while 
2NL0 correction to the ratio r is large and important. 

The rather small difference in r values results in quite noticeable dis- 
agreement of the slopes r'. Theoretical estimates can be shown IQ to be 
quite predictive for the ratio of the slopes of multiplicities but it is much less 
reliable to use the perturbative estimates even at Z°-energy for such quan- 
tities as the slope of r or the ratio of slopes of logarithms of multiplicities 
(the logarithmic slopes). Much higher energies are needed to do that. Thus 
the values of r and/or of the logarithmic slopes can be used to verify the 
structure of the perturbative expansion. 

We demonstrate it here on the example of the slope value. The slope r' 
is extremely sensitive to higher order perturbative corrections. The role of 
higher order corrections is increased here compared with r because each nth 
order term proportional to 7q gets an additional factor n in front of it when 
differentiated, the main constant term disappears and the large ratio r2/ri 
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becomes crucial: 



(11) 



where the relation 7q ?a — i?7o(l+-Bi7o) has been used with B = Po/SNc] Bi = 
(3i/4Nc(3o. The factor in front of the bracket is very small already at present 
energies: BtqTi 0.156 and 70 ~ 0.5. However, the numerical estimate of r 
is still indefinite due to the expression inside the brackets. Let us note that 
each differentiation leads to a factor as or 7q, i.e., to terms of higher order. 
For values of ri, r2, r-j, tabulated above [uf = 4) one estimates 2r2/ri ~ 4.9, 
(Sra/ri) ^ 1.5). The simplest correction proportional to 70 is more than 
twice larger 1 at energies studied and the next one is about 0.4. Therefore 
the ever higher order terms should be calculated for the perturbative values 
of r' to be trusted. The slope r' is equal to for a fixed coupling constant. 

The higher order terms are also important for the moments of the multi- 
plicity distributions [Q. The normalized second factorial moment F2 defines 
the width of the multiplicity distribution. 

The asymptotical (70 0) values of and F2 are different: 

4 7 

T?G _ Z _ ^ (-in\ 

At energy the widths of the distributions are smaller 

F^ ^1.12, F^ ^ 1.34. (13) 

but still much larger than their experimental values 1.02 and 1.08, corre- 
spondingly. The rather large difference of the perturbative (|13D and exper- 
imental values at indicates that moments of the distributions should be 
sensitive to corrections. The conclusions about the third moments are similar. 
Nonetheless, the computer solution of the QCD equation happened to 

be quite successful in fitting experimental data even for higher moments and 
their ratios Hg introduced in |]ll| . It shows that the role of conservation laws 



treated approximately in the analytical approach and accurately accounted 
in computer calculations becomes more important for higher moments. 

Thus it is shown that the analytical perturbative approach is quite suc- 
cessful in demonstrating that all features of QCD predictions about multi- 
plicities of quark and gluon jets correspond to the general trends of exper- 
imental data. Some disagreement at the level of 15-20% is understandable 
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due to incomplete account for the energy-momentum conservation in such 
an approach. Further accurate computer solutions are needed to check if 
these trends persist at the higher precision level. They would take into 
account some non-perturbative effects as well. The dipole model of QCD 
which includes implicitly the non-perturbative string effects has been devel- 
oped [0, |13[ in various versions. It describes experimental data even with 



higher precision at present energies [ p^ . 
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